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ABSTRACT 


This  report  summarizes  the  results  of  a  study  on  the  control  and 
stabilization  of  linear  and  nonlinear  distributed  systems  covering  the 
period  from  April, 1986  through  December , 1988.  The  specific  areas  of  this 
study  consist  of  (i)  stability  enhancement  of  distributed  systems 
describable  by  abstract  linear  evolution  equations,  (ii)  a  translation 
invariant  approach  to  stability  and  stabilizability,  (iii)  stabilizability 
of  bilinear  systems,  (iv)  application  of  computer  vision  in  static  shape 
estimation,  control,  and  failure  detection  in  elastic  systems,  and  (iv) 
stabilization  and  control  of  distributed  systems  with  time-dependent 
spatial  domains  with  application  to  aeroelastic  systems  with  extendible 
lifting  surfaces. 
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1.  INTRODUCTION 

The  use  of  feedback  controls  for  stabilizing  the  motions  of  distributed 
systems  has  been  investigated  in  recent  years.  Despite  the  extensive 
results  in  this  area,  it  is  still  not  completely  clear  how  to  find 
stabilizing  feedback  corttrols  for  distributed  systems,  particularly 
when  they  are  nonlinear.  In  this  study,  we  focus  our  attention  on  the 
stabilization  and  control  of  certain  classes  of  distributed  systeas 
describable  by  abstract  linear  evolution  equations,  bilinear  systeas,  and 
systems  involving  time-dependent  spatial  domains.  The  latter  type  of 
systems  arise  in  aeroelastic  systems  with  extendible  lifting  surfaces,  and 
also  in  large  space  structures  with  deployable  structural  members. 
Another  aspect  of  this  study  pertains  to  the  use  of  computer  vision  for 
shape  estimation,  control  and  failure  detection  of  deformable  bodies  such 
as  large  elastic  structures.  j)The  main  objective  of  this  study  is  to 
develop  applicable  mathematical  theories  and  at  the  Bame  time  study 
specific  systems  arising  from  realistic  aerospace  applications.  '  r 

2.  RESEARCH  SUMMARY 

This  report  covers  the  period  from  April, 1986  through  December , 1988. 
The  results  will  be  described  categorically  below: 


2.1  Stability  Enhancement  by  State  Feedback: 

Let  (A,B)  denote  the  distributed  system  described  by  the  abstract 
equation : 

x  =  Ax  +  Bu,  (1) 


where  A  is  the  generator  of  a  strongly  continuous  semigroup  T(t),  t*  0, 
i.e.,  of  the  class  C  ,  of  bounded  linear  operators  over  a  Hilbert  space  H, 

O 

and  B  is  a  bounded  linear  operator  from  a  second  Hilbert  space  V  to  H. 
Let  F:  H  -»  U,  u=Fx,  be  a  bounded  linear  feedback  so  that  (1)  becomes: 

x  =  (A  +  BF)x.  (2) 


Then,  since  BF  is  bounded  linear,  A  +  BF  also  generate;,  a  C  semigroup 

O 

S(t),  t  i  0,  (say),  over  H. 

Now,  if  the  open-loop  system  ( A , B )  is  not  stable  and  if  a  state  feedback 
F  can  be  found  so  that  the  closed-loop  system  (2)  is  stable,  then  we  say 
that  (A,B)  is  stabilizable  by  F.  This  is  the  well-known  stabilizability 
problem  in  Control  Theory.  Our  investigation  here  is  somewhat  more 
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general.  It  was  motivated  by  the  fact  that  many  physical  systems  are 
inherently  stable.  But  their  response  may  be  too  "slow".  Thus  we  wish  to 
use  state  feedbacks,  not  only  for  stabilizability ,  but  also  for  "speeding 
up"  the  response  of  a  system. 

We  began  with  the  class  of  conservative  systems,  e.g.,  those  described 
by  the  abstract  wave  equation.  For  this  class,  the  semigroup  T(t),  t  i  0, 
is  norm  preserving:  II T ( t ) x II  =  II x II  for  all  x  in  H.  Thus  a  conservative 
system  is  not  s(strongiy)-stable  on  all  the  state  space  H:  T(t)x  -»  0,  t  •* 
®,  for  each  x  in  H.  It  is  also  not  e(exponentially)-stable,  i.e.,  there 
are  no  constants  M  ^  1  and  a  >  0  so  that:  llT(t)ll  £  M  exp(-at),  t  £  0. 
Thus,  the  natural  question  to  ask  is  when  does  a  norm  preserving  semigroup 
T ( t ) ,  t  £  0,  w(weakly  (-stable:  [T(t)x,y]  -»  0,  t  -»  ®>,  for  x  and  y  in  H? 

Necessary  and  sufficient  conditions  for  a  Cq  semigroup  to  be  e-stable 
were  given  by  Datko.  They  are  indeed  extension  of  the  well-known  Lyapunov 
Theorem  for  stability  of  matrices.  However,  as  far  as  we  know,  there  are 
no  Lyapunov  type  theorems  for  weak  and  for  strong  stabilities  of  Hilbert 
space  semigroups. 

Since  a  norm  preserving  semigroup  belongs  to  the  class  of  contraction 
semigroups:  II T ( t ) II  £  1,  t  *  0,  we  began  in  [Rl]  by  obtaining  a  Lyapunov 
type  condition  for  w-stability  of  contraction  semigroups.  Namely,  for  a 
contraction  semigroup  T(t),  t  £  0,  with  generator  A  in  H,  there  always 
exists  a  non-negative  operator  P  such  that:  2Re  [PAx,x]  £  0  for  x  in  the 
domain  D( A)  of  A.  In  particular,  if  P  is  positive  then  the  semigroup  is 
weakly  stable.  Moreover,  if  the  resolvent  of  A  is  compact  —  which  is  the 
case  in  most  physical  systems  —  then  our  conditions  also  result  in  strong 
stability  of  the  semigroup. 

Returning  to  the  case  of  a  conservative  system.  The  key  idea  is  to  use 
a  feedback  so  that  the  closed-loop  semigroup  S(t),  t  i  0,  is  a  contraction 
semigroup.  One  such  feedback  is  -B*,  or  more  generally,  -Q2B  where  Q  is 
a  self-adjoint  operator  on  U.  The  latter  has  the  advantage  that  if  B  is 
not  compact,  we  can  still  have  a  compact  feedback  by  requiring  Q  to  be 
compact.  Using  our  results  stated  above,  we  were  able  to  show  that  if  A 
has  compact  resolvent,  and  if  0  is  not  in  the  spectrum  of  BB*,  then  the 
feedback  -B*  will  strongly  enhance  the  system  stability.  Stabilizability 
of  conservative  systems  was  also  studied,  using  the  same  Lyapunov  type 
results. 
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A  rather  interesting  outcome  of  our  work  is  the  fact  that  if  the 

resolvent  of  the  generator  A  of  a  conservative  system  is  compact,  then  it 
is  always  weakly  unstable  on  all  of  its  state  space.  This  explains  why, 
for  such  a  system,  controllability  is  required  on  all  the  state  space 
before  stabilizability  or  stability  enhancement  can  be  achieved.  Behavior 
of  the  modes  of  a  conservative  system  under  the  influence  of  a  dissipative 
feedback  was  also  studied.  Conditions  for  the  modes  of  the  closed-loop 
system  to  be  orthogonal  were  also  obtained. 

In  order  to  develop  Lyapunov  type  results  for  weak  and  strong 

stabilities,  we  studied  in  [R2]  the  operator  equation: 

[PAx,x]  +  [x,PAx]  =-[Rx,x]  for  x  £  D(A).  (3) 

Here  the  operator  R  is  required  only  to  be  nonnegative.  This  equation  is, 
of  course,  a  generalization  of  the  Lyapunov  equation  —  in  which  R  is 

either  the  identity  operator  or  is  strictly  positive,  i.e.,  R  <:  al  for 

some  a  >  0.  Necessary  and  sufficient  conditions  for  the  existence  of  a 
unique  non-negative  solution  P  were  obtained.  The  important  case  in  which 
R  =  BB*  was  studied.  Sufficient  conditions  for  weak  stability  of  a 
uniformly  bounded  C  semigroup  were  then  obtained.  These  conditions,  as 

O 

far  as  we  can  tell,  have  not  been  given  before.  It  is  of  interest  to  note 

•ft 

that  when  R  =  BB  our  conditions  becomes:  (i)  equation  (3)  admits  a  unique 

solution  P  2  0  if  and  only  if  the  integral  fQ  llB  T ( t ) x H 2  dt  converges  for 
every  x,  and  (ii)  if  the  pair  (A*,B)  is  controllable  then  P  is  positive 
and  the  uniformly  bounded  semigroup  is  weakly  stable.  Condition  (i)  is 
clearly  a  Lyapunov  type  condition,  while  (ii)  is  needed  for  the  positivity 
of  P.  Thus,  unlike  exponential  stability,  controllability  is  required  in 
this  case.  These  two  conditions  were  used  in  [R2].  for  stabilizability 
of  uniformly  bounded  semigroups.  Note  that  if  the  uniformly  boundedness 
condition  is  not  required,  then  the  above  conditions  imply  that  the 
semigroup  is  a  quasi-affine  transform  of  a  contraction  semigroup.  Hence 
we  can  show  that  it  is  only  weakly  stable  on  a  dense  subspace.  This  is 
precisely  the  case  when  one  uses  the  feedback  -B  P  where  P  £  0  satisfies 
the  steady-state  Riccati  equation:  [PAx,x]  +  [x,PAx]  -  [PBB  Px,x]  +  [Rx,x] 
=  0  for  x  in  the  domain  of  A,  and  for  some  R  *  0.  The  above  suggests  a 
general  relationship  between  Lyapunov  type  stability  on  the  one  hand  and 
quasi-affine  transforms  of  contraction  semigroups  on  the  other  hand. 
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Results  on  Lyapunov  type  strong  stability  were  obtained  in  [R3].  The 
approach  here  is  again  via  equation  (3).  It  turns  out  that  if  the 
solution  P  of  the  equation  is  strictly  positive,  then  the  semigroup  is 
strongly  stable.  Several  conditions  for  strict  positivity  of  P  were 
obtained.  We  must  note  that  exponential  stabilizability  is  the  most 
desirable  property,  but  this  unfortunately  is  not  always  possible!  This 
fact  emphasizes  the  need  for  feedback  controls  which  will  either  strongly 
or  weakly  stabilize  or  enhance  the  stability  of  a  system.  These  were 
studied  in  [R4].  Much  of  prior  work  in  stabilizing  controls  for 
distributed-parameter  systems  has  emphasized  exponential  stabilizability. 
Compared  with  exponential  stability,  weak  and  strong  stabilities  are  less 
desirable  properties.  However,  there  are  situations  under  which  infinite 
dimensional  systems  cannot  be  exponentially  stabilized.  In  such  cases,  we 
propose  weak  and  strong  stabilities.  One  approach  to  the  stabilization  of 
finite  dimensional  systems  and  exponential  stabilization  of  infinite 
dimensional  systems  has  been  the  use  of  Lyapunov  type  functionals.  This 
is  one  technique  which  we  developed  and  extended  here,  to  provide  new 
conditions  for  strong  and  weak  stabilities.  We  presented  a  new  functional 
in  order  to  obtain  conditions  for  a  certain  semigroup  to  be  strongly 
stable.  This  functional  also  suggested  an  inequality  relation  which,  if 
satisfied  will  guarantee  weak  stability  of  uniformly  bounded  semigroups. 
Since  stable  semigroups  are  uniformly  bounded  and  since  this  condition  is 
important  in  verifying  stability,  we  examined  this  phenomenon.  Some  new 
conditions  were  found  for  a  uniformly  bounded  semigroup  to  remain 
uniformly  bounded  under  suitable  bounded  perturbations. 

In  [R5],  strong  stability  of  distributed  systems  was  studied  by 

characterizing  strongly  stable  semigroups  of  operators  associated  with 

distributed  systems.  Main  emphasis  is  on  contractive  systems.  Three 

different  approaches  to  characterization  of  strongly  stable  contractive 

semigroups  are  developed.  The  first  one  is  an  operator  theoretical 

approach.  Using  the  theory  of  dilations,  it  was  shown  that  every  strongly 

stable  contractive  semigroup  is  related  to  the  left  shift  semigroup  on  an 
2 

L  space.  Then,  a  decomposition  for  the  state  space  which  identifies 
strongly  stable  and  unstable  states  was  introduced.  Based  on  this 
decomposition,  conditions  for  a  contractive  semigroup  to  be  strongly 
stable  were  obtained.  Finally,  extensions  of  Lyapunov  equation  for 
distributed  systems  are  investigated.  Sufficient  conditions  for  weak  and 
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strong  stabilities  of  uniformly  bounded  semigroups  were  obtained  by 
relaxing  the  equivalent  norm  conditions  on  the  right-hand-side  of  the 
Lyapunov  equation. 

These  characterizations  are  then  applied  to  the  problem  of  feedback 
stabilization.  First,  it  was  shown  via  the  state  space  decomposition  that 
under  certain  conditions  a  contractive  system  (A,B)  can  be  strongly 
stabilized  by  the  feedback  -B*.  Then,  application  of  the  extensions  of 
the  Lyapunov  equation  results  in  sufficient  conditions  for  weak,  strong 

and  exponential  stabilization  of  contractive  systems  by  the  same  feedback 

*  • 

-B  .  Finally,  it  was  shown  that  for  a  contractive  system  x  =  Ax  +  Bu, 

where  B  is  any  linear  bounded  operator,  there  is  a  related  linear 

quadratic  regulator  problem  and  a  corresponding  steady-state  Riccati 

equation  which  always  has  a  bounded  nonnegative  solution. 

2.2  A  Translation  Invariant  Approach  to  Stability  and  Stabilizability: 

This  study  was  motivated  by  the  following  observations: 

Let  A  be  the  generator  of  a  contractive  semigroup  T(y),  t  i  0,  over  a 
Hilbert  space  H.  Then  for  x  in  the  domain  D(  A)  of  A: 

[Ax,x]  +  [x,Ax]  S  0,  (4 ) 


and 

oo 

llxll?  -  lim  llT(  t  )xll 2  =  I  -2Re  [AT(  t  )x,T(  t  )x]  dt,  (5) 

t-*00  0 

where  the  limit  always  exists,  since  the  function  II T ( t ) x It  is 

nonincreasing.  Define  on  D(  A)  a  new  norm  It'll  by 


00 

xll2  =  I  II T  ( t )  x  H  2  dt  for  x  e  D(  A). 

J0 


(6) 


Let  K  be  the  completion  of  D( A)  in  the  new  norm  —  modulo  the  zero 
vectors.  Then  (6)  simply  means  that,  for  x  in  D{ A):  T(t)x  belongs  to  the 


space  L  ([0,®);K).  Define  an  operator  V:  D{ A) 
(Vx)(t)  =  T(t)x  for  t  i  0. 


L  ( [0,®) ;K)  by 


Then  it  is  plain  from  (6)  that  V  is  an  isometry  on  D{  A).  Therefore,  since 
D( A)  is  dense  in  H,  we  can  extend  V  by  continuity  to  all  of  H.  More  is 
true.  We  have,  for  x  in  H: 

VT ( t ) x  =  T( t )T{ t ) x  =  T ( t+t )x  for  t,T  i  0. 
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Therefore 


VT(t)x  =  L(T)T(t)x  =  L(T)x  for  x  in  H,  (7) 

where  L(t),  t  2:  0,  is  the  Left  Translation  semigroup  over  the  space 

2 

L  ( [ 0 ,<*>) ; K ) .  We  therefore  conclude  that  the  closed  subspace  VH  is 
invariant  for  the  semigroup  L(t),  t  i  0,  and  the  strongly  stable 
contraction  semigroup  T(t),  t  £  0,  to  its  invariant  subspace  VH.  It  is 
clear  that  if  a  semigroup  is  unitarily  equivalent  to  a  part,  i.e.,  the 
restriction  to  an  invariant  subspace,  of  the  Left  Translation  semigroup, 
then  it  is  strongly  stable. 

We  have  in  the  above  obtained  a  Translation  Invariant  Representation  for 
a  strongly  stable  contraction  semigroup. 

In  [R6],  we  were  able  to  obtain  a  general  Translation  Invariant 
Representation  for  a  contraction  semigroup  by  using  the  fact  that  the 

contractions  T(t)*T(t),  t£  0,  converge  in  the  strong  operator  topology  to 

2 

a  positive  contraction  C  (say).  Therefore,  from  (5),  for  x  in  D( A): 

a 

llxll2  -  It Cx II 2  =  I  II T ( t ) x II 2  dt.  (8) 

J0 

Therefore,  for  x  in  D( A): 

oo 

llxll?  =  I  llT ( t ) x II 2  dt  +  II Cx II 2 . 

J  0 

From  this  we  can  show  that  T(t),  t  £  0,  is  unitarily  equivalent  to  the 

restriction  of  the  semigroup  L(t)  ®  W(t),  t  i  0,  to  its  invariant  subspace 

2  — 

VH,  where  V  is  now  the  isometry  from  H  to  the  space  L  ([0,®);K)  ®  CH 

defined  by  (Vx)(t)  =  T(t)x  ®  W(t)Cx,  and  W(t),  t  it  0,  is  an  isometric 

semigroup  over  CH  defined  by  W(t)Cx  =  CT(t)x.  What  interesting  is  that  in 

the  process  of  deriving  the  Translation  Invariant  Representation  we  obtain 

yet,  another  Lyapunov  type  result  for  w-stability  of  contraction 

2 

semigroups.  This  results  directly  from  (8)  by  defining  P  =  I  -  C  .  Then 
(8)  becomes 

<!• 

[Px,x]  =  [  II T ( t ) x II 2  dt  for  x  e  Z>( A).  (9) 

J0 

From  which  it  is  easy  to  see  that  the  semigroup  is  w-stable  as  soon  as  P 
is  positive. 

The  Translation  Invariant  Representation  is  not  only  limited  to  the 
class  of  contraction  semigroups.  We  were  able  to  show  that  a  class  of 
exponentially  stable  C  semigroups  also  admits  such  a  representation.  To 
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see  this  we  only  have  to  note  that  if  T(t),  t  £  0,  is  exponentially 
stable,  then  there  exists  a  positive  operator  P  such  that,  for  every  x: 

CD 

II T ( t ) x II 2  dt.  (10) 

0 

2 

This,  as  in  the  above,  shows  that  T(t)x  belongs  to  the  space  L  {[0,®);H). 

2 

Now  define  a  new  norm  11*11  by:  II x II p  =  [Px,x],  and  suppose  that  it  is 

equivalent  to  the  original  norm  11*11.  Then  it  is  plain  from  (10)  and  from 

the  above  discussion  that  the  semigroup  is  unitarily  equivalent  to  the 

2 

left  shift  semigroup  over  the  space  L  ([0,®);H).  Conditions  for  the  two 
norms  II -  II p  and  11*11  to  be  equivalent  have  been  given  by  Pazy.  It  turns  out 
that  they  are  equivalent  for  the  "hyperbolic"  case,  but  not  for  the 
"parabolic"  one.  We  must  note  that  if  the  generator  A  is  bounded  and  the 
semigroup  exp(At),  t  ^  0,  is  exponentially  stable,  then  the  two  norms  are 
equivalent.  This  immediately  shows  that  stability,  stabilizability  and 
stability  enhancement  in  finite  dimensional  case  can  all  be  studied  from 
the  Translation  Invariant  Representation  viewpoint. 

2.3  Stabilizability  of  Bilinear  Systems: 

Consider  the  bilinear  system 

x(t)  =  Ax  ( t. )  +  B(x(t),u(t)),  t  £  0, 

where  A  is  the  generator  of  a  linear  Cq  semigroup  over  a  Hilbert  space  H, 
and  B ( •  ,  -  )  is  a  bilinear  form  from  the  product  space  H  *  U  to  H.  Our 
problem  here  is  to  stabilize  such  a  system  by  state  feedback. 

In  [R7],  we  studied  the  class  of  "dissipative"  bilinear  systems,  i.e., 
those  for  which  the  gener<ifor  A  is  dissipative.  The  idea  here  is  to  use  a 
state  feedback  so  that  the  closed-loop  system  is  also  dissipative.  We 
studied  the  case  in  which  B(*,*)  is  of  the  form  B(x(t),u(t))  =  u(t)Qx(t), 
where  Q  is  a  bounded  linear  operator  on  H.  A  suitable  feedback  was 
obtained,  and  the  system  is  weakly  stabilizable  as  soon  as  it  is  "weakly 
controllable",  i.e.,  [ B*T{ t )x ,T( t )x ]  =0,  t  £  0,  •*  x  =  0. 

Next,  we  consider  the  class  of  Semilinear  Dissipative  Systems: 

x(t)  =  Ax(t)  +  B(x(t))v(t). 

Sufficient  conditions  for  this  class  of  systems  to  be  strongly 
stabilizable  by  the  feedback  -B*  were  obtained.  Namely, 

(x:  B*(T(t)x)T(t)x  =  0,  t  i  0|  =  (0). 
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This  is  the  exact  analog  of  the  weak  controllability"  condition  above. 
However,  in  this  case,  stronger  results  was  obtained.  This  requires  the 
knowledge  of  the  "isometric"  subspace  of  the  contraction  semigroup 
generated  by  A,  which  is  also  the  case  for  the  class  of  linear  dissipative 
systems. 

2.4  Application  of  Computer  Vision  in  Static  Shape  Estimation,  Control 
and  Failure  Detection  in  Elastic  Systems: 

This  study  was  motivated  by  the  observation  that  a  computer-vision 
system  is  capable  of  captur.ng  the  instantaneous  shape  of  a  deformable 
solid  body  over  its  entire  spatial  domain.  Moreover,  the  sole  dependence 
of  vision  data  on  the  local  surface  properties  of  the  body  eliminates  the 
need  of  a  precise  mathematical  model  of  the  system  under  consideration. 
These  features  make  the  vision  system  particularly  attractive  for  the 
feedback  control  of  elastic  systems  involving  deformable  solid  surfaces. 
Moreover,  computer-vision  is  capable  of  detecting  large  deformations  and 
abnormalities  in  an  elastic  structure.  This  capability  may  be  used  to 
detect  failures  and  abnormal  conditions  in  elastic  structures.  The  main 
objective  of  this  study  is  to  develop  applicable  mathematical  theories  for 
shape  estimation,  control  and  failure  detection  of  elastic  systems  arising 
in  aerospace  applications. 

In  order  to  gain  an  understanding  of  the  basic  problems  and  issues  in 
the  application  of  computer-vision  systems  for  shape  estimation,  we  began 
in  1986  with  an  experimental  study  on  static  shape  estimation  using  a 
small  scale  model  of  a  spaceborne  antenna  consisting  of  two  flexible  ribs 
interconnected  by  elastic  strings.  A  crude  CCD  camera  with  binary  grey 
levels  was  used  in  the  experiment.  This  study  led  to  the  formulation  of 
physically  meaningful  mathematical  problems  in  shape  estimation.  The 
preliminary  results  of  this  study  were  presented  at  the  1986  IFAC 
Symposium  on  the  Control  of  Distributed-Parameter  Systems  [ R8  ] . 

During  the  period  1986-1987,  with  the  collaboration  of  the  Jet 
Propulsion  Laboratory  (JPL)  in  Pasadena,  we  studied  the  application  of 
computer  vision  in  static  shape  estimation  for  elastic  systems  with  the 
aid  of  the  a  reduced  scale  model  of  a  large  spaceborne  antenna  at  the  JPL 
Ground  Antenna  Test  Facility.  This  model  consists  of  twelve  flexible  ribs 
attached  to  a  rigid  hub  which  is  connected  to  another  rigid  ring  through  a 
gimbal  system.  The  entire  antenna  model  has  two  degrees  of  freedom 
associated  with  the  rigid  body.  The  ribs,  hub  and  the  rings  were  coated 
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with  light  reflecting  paint  so  they  could  be  viewed  by  Beans  of  a  computer 
vision  system.  Our  initial  objective  was  to  estiaate  the  displacement  of 
each  rib  relative  to  a  suitable  hub  coordinate  system. 

The  shape  estimation  problem  can  be  decomposed  into  three  subproblems: 
(PI)  estimation  of  hub  attitude  relative  to  a  specified  inertial  frame; 
(P2)  estimation  of  rib  displacement  relative  to  the  hub  coordinate  system, 
and  (P3)  reconstruction  of  the  reflector  surface  from  the  estimated  rib 
shape. 

Each  of  the  foregoing  problems  can  be  considered  as  a  special  case  of 
the  following  general  mathematical  problem: 

Let  T  be  a  nonempty  compact  subset  of  the  3-dimensional  real  coordinate 

3 

space  R  representing  the  deformable,  stationary  solid  object  under 
observation  (eg.  antenna  hub,  rib  or  reflector).  The  boundary  of  T  is 
denoted  by  dF.  To  view  dF  by  means  of  a  computer  vision  system,  we  cover 
dr  with  a  light-reflecting  coating  or  introduce  appropriate 
light-reflecting  markings  over  a  subset  dr  of  dF.  We  assume  that  dr  or 

M 

its  markings  are  illuminated  by  an  appropriate  light  source.  The 

reflected  light  is  collected  by  the  lens  system  of  the  camera  and 
projected  into  an  image  plane.  This  image  formation  process  can  be 
described  in  abstract  mathematical  terms  as  follows: 

Let  5  be  a  family  of  planes  11.  in  R  .  The  camera’s  image  plane 
corresponding  to  a  particular  camera  orientation  is  represented  by  a 

compact  subset  of  D.  for  some  i.  The  domain  of  the  image  of  dr  (resp. 

dr  ), denoted  by  (dT)  (resp.  (dT  ))  is  specified  by  P .  ( dr )  n  y.  (resp. 

n  1  1  H  1  1 

3 

P  (dr  )  n  y  ),  where  P  represents  the  image  projection  operator  from  R 
— *  n,.  Note  that  for  perspective  image  projection,  P.  is  nonlinear. 
However,  by  using  the  "homogeneous  coordinates",  P.  can  be  represented  by 
a  linear  transformation  P(  defined  on  the  augmented  4-dimensional  space 
R4.  The  light  intensity  of  the  image  V.{dT)  (resp.  9 .  (dT  ))  can  be 

i  in 

described  by  a  real-valued  nonnegative  function  f.  defined  on  5. (dr) 

(resp.  ST, (dr  )).  The  form  of  f.  depends  on  the  light  reflecting 

I  n  i 

properties  of  dr  (or  its  markings  dT  ) ,  and  the  viewing  geometry. 

n 

In  setting  up  the  computer  vision  system  for  determining  the  shape  of 
the  object  T,  it  is  of  importar~e  to  choose  the  camera  position  and 
orientation  such  that  the  object’s  surface  dr  can  be  determined  from  its 
image  f!  (dT).  This  problem  can  be  stated  as  follows: 

(PI)  Camera  Positioning  Problem:  Given  f  and  T,  find  a  plane  fl.e  5  such 
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that  ar  can  be  determined  completely  from  its  image  ?.(dr).  In  the  case 
of  a  dual-camera  system,  we  wish  to  find  a  pair  (II., 11^),  i  t  j,  11^,  Il^e  9, 
such  that  ar  can  be  determined  completely  from  the  corresponding  image 
pair  (^.(ar),s^(a r ) ) . 

The  foregoing  problem  is  a  pure!’  geometric  one.  Suppose  that  this 

problem  has  been  solved.  We  can  now  state  a  shape  estimation  problem 

under  the  assumption  that  a  perfect  image  of  the  object  is  available.  In 

what  follows,  we  shall  drop  the  subscript  i  with  the  understanding  that 

suitable  camera  position  and  orientation  have  been  selected. 

(P2)  Idealized  Shape  Estimation  Problem  (Single  Camera  Vision  System): 

Given  an  image  intensity  function  f  defined  on  9(dT)  corresponding  to  a 

~  3 

specified  camera  orientation , f ind  a  set  3r  c  R  which  best  approximates  dr 
in  some  sense. 

In  an  actual  computer  vision  system,  the  actual  image  of  dr  is  formed 

through  both  spatial  discretization  and  intensity  quantization  processes. 

Let  the  image  plane  9  be  partitioned  into  compact  subsets  A^  (pixels)  such 

that  9  =  u  A  .  Let  the  characteristic  function  of  A  be  denoted  by  S  . 

i  J  j  -  J 

Then  the  spatially  discretized  image  intensity  function  is  given  by  f  = 
E.fS^  defined  on  J .  Let  LN=  , . . .  )  c  R  denote  the  finite  set  of 

quantized  light  intensity  levels.  We  introduce  an  intensity  quantization 
rule  specified  by  the  function  Q  which  maps  each  f  =  fS .  into  L  .  Thus, 

j  j  N 

the  quantized  image  intensity  corresponding  to  pixel  A^  is  given  by  f^  = 
QfS^.  Now,  a  discretized  version  of  the  shape  estimation  problem  can  be 
stated  as  follows: 

(P3)  Discrete  Shape  Estimation  Problem  :  Given  a  spatially  discretized 

intensity-quantized  image  of  dT  represented  by  the  intensity  sequence 

3 

{ f ^ } ,  find  a  set  dr  c  r  which  best  approximates  dr  in  some  sense. 

In  physical  situations,  the  image  of  the  object  of  interest  may  be 
corrupted  by  noise  induced  by  spatial  discretization  and  intensity 
quantization,  and  by  undesirable  reflected  light.  Moreover,  the  total 
image  may  also  contain  images  of  background  objects.  For  example,  in 
estimating  the  rib  shape  of  the  JPL  antenna  model  by  means  of  computer 
vision,  the  hub  and  inner  ring  may  also  appear  in  the  image.  It  is 
necessary  to  extract  the  rib-image  data  from  the  total  image  so  that  only 
relevant  data  are  used  in  rib-shape  estimation.  To  reduce  the  effect  of 
noise,  some  form  of  image  preprocessing  may  be  introduced.  The  extraction 
of  the  image  of  the  object  of  interest  from  the  total  image  involves  image 
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segmentation  and  object  recognition.  In  the  JPL  antenna  experiment,  this 
problem  is  solved  by  introducing  appropriate  localized  illumination  of  the 
object. 

During  the  period  covered  by  this  report,  explicit  forms  for  various 
mappings  (i.e.  P  ,  f.  etc.)  associated  with  problems  (P1)-(P3)  were 
obtained.  Algorithms  for  shape  estimation  based  on  least-square  criteria 
were  derived  and  implemented.  These  algorithms  were  tested  using  real 
data  obtained  from  an  experiment  performed  at  the  JPL  Ground  Antenna  Test 
Facility  during  the  summer  of  1987.  The  results  were  described  in  a  UCLA 
Engineering  report  [RIO]. 

After  completion  of  the  aforementioned  experimental  work,  we  extended 
our  study  to  the  detection  of  failures  and  abnormal  conditions  in  large 
spaceborne  antenna  systems.  The  basic  idea  is  as  follows: 

N  N  A  A 

Let  dr  and  dr  (resp.  dr  and  dr  )  denote  respectively  the  boundary 

n  H 

surface  of  the  light-reflecting  markings  of  the  spatial  domain  of  the 
antenna  system  in  its  normal  state  (resp.  actual  state).  The  images  of 
dr*  and  df^ ,  i=  N  or  A,  under  the  projection  operator  P  in  the  camera’s 

image  plane  9"  are  denoted  by  S^dT1)  =  P(dr‘)  n  JP  and  ^(dT*  )  =  P(dTl  )  n  ? 

M  M 

respectively.  The  projection  operator  P  depends  on  the  camera  parameters 
and  orientation  which  are  assumed  to  be  fixed  at  all  times.  To  fix  ideas, 
we  assume  that  the  camera  has  binary  intensity  levels.  Thus,  the  images 
2f(drJ*)  and  E7(dT*)  are  compact  subsets  of  ?  c  R2. 

M  n 

Now,  the  basic  approach  for  detection  of  failure  and  abnormal  condition 

is  as  .ollows:  First,  a  set  of  reference  images  9  (dr  ) ,  jrl , , , ,  ,K, 

J  H 

corresponding  to  the  normal  state  of  the  antenna  structure  are  obtained  by 
means  of  a  computer  vision  system.  Then,  the  state  of  the  antenna 
structure  is  monitored  periodically  by  the  same  computer  vision  system. 
The  resulting  images  are  compared  with  their  corresponding  reference 
images  and  their  deviations  are  measured  in  terms  of  a  suitable  index. 

N  A 

Here,  the  deviation  between  two  images  £7.  (dr  u)  and  (7.(dr  )  is  taken  to  be 

j  M  j  M 

the  Hausdorff  metric  defined  by 

p(fJ  .(dr"),sr  (dr*))  =  max{max[p(r,£7.(ar||)),max[p(r’,£7  (dr*))}}, 

jnjn  jn  jn 

where  p(r,A)  =  infjllr  -  r’ll,  r’  e  a).  The  antenna  is  said  to  be  in  the 
"failed"  or  "abnormal"  state  when  the  deviation  between  the  reference  and 
actual  images  exceeds  a  predetermined  threshold  level.  We  may  also 
introduce  various  weights  in  the  Hausdorff  metric  to  emphasize  the 
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relative  importance  of  deformations  in  different  parts  of  the  antenna 
structure.  Alternatively,  the  same  approach  may  be  used  to  monitor  the 
state  of  certain  critical  portions  of  the  antenna  by  localizing  the  camera 
image  field  to  these  portions  of  the  antenna.  Thus,  the  nature  of  the 
failure  can  be  diagnosed.  An  obvious  advantage  of  the  foregoing  approach 
is  that  the  Hausdorff  metric  is  capable  of  taking  into  account  failures 
due  to  breakage  of  the  structural  members  and  punctures  of  the  reflector 
surface  due  to  meteorite  impact.  However,  for  a  given  antenna  structure, 
there  may  exist  structural  deformations  and  movements  which  leave  the 
image  invariant.  This  difficulty  may  be  remedied  partially  by  using  a 
multi-camera  system.  In  physical  situations,  it  is  of  importance  to 
determine  the  sensitivity  and  robustness  of  the  proposed  failure  detection 
scheme  with  respect  to  variations  in  threshold  and  illumination  levels, 
and  camera  parameters.  Mathematically,  this  task  requires  the  computation 
of  the  Gateaux  differentials  of  the  Hausdorff  metric  with  respect  to 
certain  parameters  associated  with  the  projection  operator  P. 

The  results  of  this  study  were  presented  at  the  Fifth  IFAC  Symposium  on 
Control  of  Distributed  Parameter  Systems  held  at  Perpignan,  France  in 
June,  1989.  [ R 1 2 ] .  These  results  also  constitute  a  part  of  the  Ph.D. 

dissertation  of  Nai-jian  Wang  [HI 3]. 

2.5.  Stabilization  and  Control  of  Distributed  Systems  with  Time-Dependent 
Spatial  Domains: 

Distributed  system  with  time-dependent  spatial  domains  arise  naturally 
in  many  physical  situations.  Often,  it  is  desirable  to  control  the 
dynamic  behavior  of  such  systems  by  varying  their  spatial  domains.  For 
example,  in  large  space  structures  with  deployable  elastic  components,  it 
is  of  interest  to  control  the  deployment  of  these  components  so  that  the 
induced  vibrations  will  be  below  a  specified  level. 

In  this  study,  the  problem  of  stabilization  and  control  of  distributed 
systems  with  time-dependent  spatial  domains  was  considered.  Here,  the 
evolution  of  the  spatial  domain  with  time  is  described  by  a  finite 
dimensional  system  of  ordinary  differential  equations,  while  the 
distributed  systems  are  described  by  first  and  second-order  linear 
evolution  equations  defined  on  appropriate  Hilbert  spaces.  The  controls 
affect  the  system  domain  motion  through  the  ordinary  differential 
equation. 

First,  results  pertaining  to  the  existence  and  uniqueness  of  solutions 
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of  the  system  equations  were  obtained.  For  the  case  with  first-order 
evolution  equations,  the  existence  result  was  proved  in  a  constructive  way 
using  finite  dimensional  Galerkin  approximations  with  an  appropriate 
time-dependent  basis.  This  approach  is  useful  in  numerical  computations. 
For  the  case  of  second-order  evolution  equations,  we  did  not  succeed  to 
prove  the  existence  of  solutions  using  the  same  approach. 

Next,  we  considered  various  optimal  control  problems  associated  with 
this  class  of  systems.  Sufficient  conditions  in  the  form  of  variational 
inequalities  for  which  the  optimal  controls  must  satisfy  were  established. 
We  also  considered  the  feedback  stabilization  of  such  systems. 
Stabilizing  control  laws  based  on  energy  inequalities  were  derived.  The 
application  of  the  results  were  illustrated  by  simple  examples. 

The  results  of  this  study  were  described  in  a  technical  report  £ R— 11], 
which  has  been  accepted  for  publication  in  the  Journal  of  Optimization 
Theory  and  Applications. 

2.6  Stabilization  and  Control  of  Aeroelastic  Systems  with  Extendible 
Lifting  Surfaces: 

This  study  was  motivated  from  the  problems  associated  with 
trans-atmospheric  vehicles  which  are  capable  of  multi-regime  flights  (i.e. 
subsonic, supersonic  and  hypersonic  regimes).  In  these  vehicles,  movable 
lifting  and  control  surfaces  may  be  introduced  to  provide  the  desired 
lifting  and  maneuverability  characteristics.  The  extension  or  contraction 
of  these  surfaces  during  flight  could  induce  undesirable  aeroelastic 
vibrations  which  must  be  stabilized  so  that  their  effectiveness  is  not 
impaired.  In  this  study,  we  considered  an  aeroelastic  cantilever  wing 
which  may  extend  or  contract  during  flight. 

The  wing  was  modelled  by  a  coupled  system  of  partial  differential 
equations  defined  on  a  time-dependent  spatial  domain  describing  the 
bending  and  torsional  motions  of  the  variable  length  wing.  Attention  was 
focused  on  the  following  questions: 

(a)  How  does  the  extending  or  contracting  motion  of  the  wing  affect  the 
aeroelastic  vibrations? 

(b)  I s  it  possible  to  stabilize  the  aeroelastic  vibrations  by 
introducing  appropriate  feedback  controls? 

(c)  How  do  we  introduce  suitable  finite  dimensional  approximations  for 
the  numerical  solution  of  the  system  equations? 

It  was  shown  that  when  the  wing  extends  with  constant  velocity,  the 
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vibrational  energy  increases  monotonically  with  time.  This  implies  that 
the  extending  motion  has  a  destabilizing  effect  on  the  aeroelastic 
vibrations.  On  the  other  hand,  when  the  wing  contracts  with  uniform 
velocity,  the  vibrational  energy  decreases  with  time  implying  a 
stabilizing  effect.  Although  a  wing  cannot  extend  indefinitely  during 
flight,  it  is  of  importance  that  the  amplitude  of  the  aeroelastic 
vibrations  induced  by  wing  extension  does  not  build  up  to  a  high  level 
during  the  extension  period.  Therefore  it  is  of  interest  to  stabilize  the 
motion-induced  instability  by  introducing  appropriate  feedback  controls. 
It  was  shown  by  introducing  a  simple  modification  of  the 
proportional-plus-rate  feedback  control  law  depending  on  the  wing-tip 
vibration  velocity  for  controlling  the  extending  motion  of  the  wing,  the 
motion-induced  instability  can  be  stabilized.  The  resulting  feedback 
system  is  in  the  form  of  a  partial  differential  equation  with  a  free 
boundary.  It  was  also  shown  that  a  useful  finite  dimensional  Galerkin 
approximation  can  be  obtained  by  introducing  a  suitable  time-dependent 
basis  corresponding  to  the  eigenfunctions  of  the  biharmonic  and  elliptic 
operators  associated  with  the  wing  at  any  fixed  time. 

The  results  of  this  study  were  presented  at  the  First  International 
Conference  in  Industrial  Applied  Mathematics  [PI].  The  details  were 
first  described  in  a  technical  report  [ R9 ]  which  was  published  later  in 
the  journal  Dynamics  and  Stability  of  Systems  in  1988. 
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